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In this paper we discuss the electronic kinetic energy in 
the ground state of the Holstein Hamiltonian in one space 
dimension using the Global-Local variational method together 
with perturbation theory at weak and strong coupling. The 
electronic kinetic energy is detailed over a broad region of 
the polaron parameter space, and is analyzed to accurately 
determine the location of the self-trapping transition. The 
self-trapping line so determined separates the small polaron 
regime from the large polaron regime, constituting a polaron 
phase diagram. 



PACS numbers: 71.38. +i, 71.15.-m, 71.35.Aa, 72.90. +y 



I. INTRODUCTION 

Despite the vast amount of effort and creativity that 
has been applied to the polaron problem for more than 
half a century, there remain significant aspects of po- 
laron structure and behavior that have defied satisfying 
explanation and quantitative description. Central among 
these is the so-called self-trapping transition. The self- 
trapping transition describes the change in polaron struc- 
ture from being large-polaron-like over one region of the 
polaron parameter space to being small-polaron-like in 
another. Perturbation theories exist providing asymp- 
totically accurate descriptions of polaron structure and 
properties on each side of this transition; however, the 
intermediate-coupling region in which the transition lies 
has proven very resistant to analysis. 

Our starting point is the 1-D Holstein Hamiltonian 

m 

H = -J^(4(a„ + i + a n -x) + huj^blfin 

n n 

-gftio 4a«0n + M , (1) 

n 

in which a\ creates a single electronic excitation in the 
rigid-lattice Wannier state at site n, and creates a 



quantum of vibrational energy Tiuj in the Einstein oscil- 
lator at site n. The hopping matrix element connecting 
nearest neighbors is given by J, and g is the electron- 
phonon coupling strength. 

The last few years have seen considerable progress in 
the application of polaron-friendly theoretical techniques 
to the Holstein model; these include variational tech- 
niques l^), cluster diagonalization j7|-^0[ dynamical 
mean field theory 0Q, perturbation theory [ fl3| [Tof . 
density matrix renormalization group (DMRG) E7ll8[ 
and quantum Monte Carlo simulations (QMC) [[19|-23f. 
While each methodology has its own merits, each also has 
its weaknesses and only few of them possess the scope to 
be accurate over broad regions of the parameter space, 
especially at the intermediate elcctron-phonon coupling 
well away from all limits. 

We approach the problem variationally using the 
Global-Local method, introduced by Brown et al. Q for 
the determination of polaron Bloch states and band en- 
ergies across the entire Brillouin zone (k £ {— tt, +tt}). 
In this paper, we shall focus on the polaron ground state 
(k = 0), and in particular on the kinetic energy in the 
ground state. Details regarding the variational method 
and its numerical implementation can be found in ref [jj) . 
The numerical data presented herein were computed on a 
I-D lattice of 32 sites with periodic boundary conditions. 



II. THE KINETIC ENERGY 

We determine the kinetic energy variationally, and 
compare our results with both weak-coupling perturba- 
tion theory and strong-coupling perturbation theory, as 
well as selected results of other methods in the intermedi- 
ate regime. In each case, the kinetic energy is computed 
as the expectation value 



E, 



a' (a 



n+l 



(2) 



in the appropriate ground state \tp). 

Our numerical results are summarized in Figure [I]. 
Each curve in Figure |l| contains 80-200 data points, with 
a higher density of points clustered in the steeper por- 
tions of each curve where the self-trapping transition is 
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expected to be found. No smoothing has been performed; 
each " curve" is a polygonal arc connecting computed en- 
ergies. Over most of the parameter space we have in- 
vestigated, computational errors are smaller than can be 
meaningfully conveyed in with any resolvable symbol; the 
principal exception is at weak coupling and smaller J 
values where decreasing sensitivity of the ground state 
energy to certain details of polaron structure eventually 
hampers convergence. 

With increasing adiabaticity, here beginning at 
J/hui ss 7, the ability of the variational method to repre- 
sent the complexity of polaron structure in the immediate 
vicinity of the self-trapping transition eventually is over- 
taxed, and discontinuities appear in estimated quantities 
such as the kinetic energy. Although the value of the ki- 
netic energy in the immediate vicinity of such anomalies 
is necessarily distorted, the location of the discontinu- 
ities continues to provide reasonable estimates for the 
location of the self-trapping transition, and outside of a 
narrow region, quantitative accuracy remains excellent. 
|0j For such reasons, we retain in the following data for 
the cases J/Tiuj = 7 and 9. 




FIG. 1. The polaron kinetic energy Ekin as a func- 
tion of the electron- phonon coupling. Curves left to right: 
J/fwj = 0.25, 0.5, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0 and 9.0. 

The weak coupling regime can be understood through 
comparison with weak-coupling perturbation theory, 
which to second order in the coupling constant yields 
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(3) 



Characteristically, the kinetic energy is a weak function 
of the electron-phonon coupling below the self-trapping 
transition, and grows increasingly weak with increas- 
ing adiabaticity. Owing to the minimal involvement 



of phonons in the polaron in this regime, the quasi- 
particle can be fairly characterized as a quasi-free elec- 
tron with a slightly reduced bandwidth. It is charac- 
teristic of weak-coupling perturbation theory, however, 
that the integration of phonons into polaron structure is 
Mnderrcpresentcd, causing the electron to appear more 
"free" than it actually is. 

The strong coupling regime can be understood through 
comparison with strong-coupling perturbation theory, 
which to second order in J (following the Lang-Firsov 
transformation p4j) yields |lqj2q] 
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J(2g 2 ) + f(g 2 )) 



(4) 
(5) 



wherein f(y) = Ei(y) — 7 — ln(y), and Ei(y) is the expo- 
nential integral. At very strong coupling, the kinetic en- 
ergy decays to zero, suggesting that the dressed electron 
becomes essentially immobile relative to the quasi-free 
electron; significant, however, is the fact that this decay 
is not ultimately exponential in the coupling constant 
as results from the small polaron approximation (first 
order Lang-Firsov), but a much weaker inverse power as 
suggested by perturbative corrections (second-order Lang 
Firsov) f2(|. 

A comparison of our numerical kinetic energies with 
weak (||) and strong ([|) coupling perturbation theories is 
shown in Figure 2. 
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FIG. 2. Comparison of Global-Local kinetic energy, E^, 
with weak-coupling perturbation theory, E^^ according to 
([]) (upper panel), and with the Lang-Firsov approximation, 
E^in according to (^) (lower panel), as functions of the elec- 
tron-phonon coupling. 

For weak electron-phonon coupling, the agreement be- 
tween our variational and perturbative kinetic energy is 
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very good up to g values quite close to the self-trapping 
transition. For this weak coupling case, perturbation the- 
ory overestimates the kinetic energy value (i.e. phonon 
involvement in the polaron is greater than is captured by 
perturbation theory). 

For strong coupling, our variational calculation ap- 
proaches the second-order Lang-Firsov result above the 
self-trapping transition; however, deviations persist sig- 
nificantly into the strong coupling regime, with larger J's 
converging to the Lang-Firsov result more slowly than 
smaller J's. Strong coupling perturbation theory system- 
atically underestimates the kinetic energy (i.e. phonon 
involvement in the polaron is weaker than is embodied in 
Lang-Firsov). 
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FIG. 3. Comparison of the kinetic energy Ekin as deter- 
mined by the Global-Local method (solid line) with weak cou- 
pling perturbation theory (dashed line), strong coupling per- 
turbation theory (dotted line), quantum Monte Carlo simula- 
tion (scatter plot), and density matrix renormalization group 
(diamonds) for J/Tilu = 1. The (+) symbol indicates the 
location of the self-trapping transition as determined using 
([]) on the GL data. QMC data consists of two statisti- 
cally-independent scans at the inverse temperature /3 = 20.0; 
data kindly provided by H. De Raedt |23|] . DMRG data 
kindly provided by E. Jeckelmann jTJ| . 



III. THE SELF-TRAPPING LINE 



diagram, and is intimately involved with, if not always 
ultimately responsible for, many of the difficulties en- 
countered in polaron theory. The physical transition is 
smooth at finite g and J/huj |27|]; however, it is common 
for approximate descriptions of the phenomenon cither 
to miss the transition completely, or to break down in 
some respect in its vicinity. 

Here, we will use the kinetic energy to define a transi- 
tion on the physical rationale that this energy should per- 
haps most directly reflect the heavily-dressed and essen- 
tially immobile character of the self-trapped state 
As our objective criterion for locating this transition, 
we associate the transition with a particular coupling 
strength, gsT(J/^), where the kinetic energy changes 
most rapidly with respect to g at fixed J/Tilo; this point 
is identified by the extremum 



d 2 Eki 



dg 2 







(6) 
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Since the sigmoidal shape of the kinetic energy is well 
resolved for all hopping integral values used in this work, 
a clear definition of the transition can be obtained from 
this criterion. A summary of our calculation is shown in 
Figure |I| 
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FIG. 4. Polaron phase diagram, showing the location of 
the self trapping transition. Dots (•) corresponds to our nu- 
merical calculation , dash ed line corresponds to the empirical 
curve gsr = 1 + y/ J/hu. 



The self-trapping transition is the more-or-less rapid 
change in polaron structure from that typical of large 
polarons (below the transition) to that typical of small 
polarons (above) as g or J are varied and is typically 
resolved as a feature that grows increasingly sharp with 
increasing adiabaticity. The self-trapping transition is 
certainly the most exotic feature of the polaron phase 



The function 

g ST = 1 + y/J/hw (7) 

describes a strictly empirical curve that we have intro- 
duced to simply, conveniently, and accurately character- 
ize the self-trapping line over the investigated regime. No 
fit has been performed to obtain this curve. 
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The accuracy with which gsr describes the appar- 
ent self-trapping line far exceeds any other estimation 
currently available. For example, one frequently finds 
the self-trapping transition associated with the apparent 
breakdown point of strong-coupling perturbation theory 
(at A = g 2 hu>/2J rj i) or with the crossing points of weak 
and strong-coupling perturbation theory (see, e.g., Fig- 
ure ||; neither of which compares as favorably as (0) over 
the investigated regime. There is reason to expect (0) to 
continue to accurately locate the self-trapping transition 
beyond the investigated regime, since it is known from 
strong-coupling perturbation theory that the transition 
is associated with the limit [EoTl 



[12] S. Chichi, F. D. Pasquale, and D. Feinberg, Europh. Lett. 
30, 151 (1995). 



lim 

g,J— >c» 



\ ST = \ c = - 



(8) 



which is satisfied by (0). 

Although the numerical data in Figure ^ were gen- 
erated from an analysis of one particular quantity, the 
kinetic energy, the same curve gsT can be expected to 
describe the appearance of the self-trapping transition in 
other physical quantities, verifying (]?]) to be a robust as 
well as accurate locator of the self-trapping transition, as 
will be demonstrated elsewhere [E8fi . 
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